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i^ ■ ABSTRACT 

C^ I The number of multi-planet systems known to be orbiting their host stars with orbital periods 

that place them in mean motion resonances is growing. For the most part, these systems are 

^SJ \ in first-order resonances and dynamical studies have focused their efforts towards understanding 

^ . the origin and evolution of such dynamically resonant commensurabilities. We report here the 

\l \ discovery of two super-Earths that are close to a second-order dynamical resonance, orbiting the 

^ ' metal-poor ([Fe/H]=-0.43 dex) and inactive G2V star HD41248. We analysed 62 HARPS archival 

r — I radial velocities for this star, that until now, had exhibited no evidence for planetary companions. 

Using our new Bayesian Doppler signal detection algorithm, we find two significant signals in the 
data, with periods of 18.357 days and 25.648 days, indicating they could be part of a 7:5 second- 

r^ . order mean motion resonance. Both semi-amplitudes are below Sms^^and the minimum masses 

of the pair are 12.3 and 8.6 M®, respectively. Our simulations found that apsidal alignment 
stabilizes the system, and even though libration of the resonant angles was not seen, the system 
is affected by the presence of the resonance and could yet occupy the 7:5 commensurability, which 

^^ . would be the first planetary configuration in such a dynamical resonance. Given the multitude 

3 ' of low-mass multiplanet systems that will be discovered in the coming years, we expect more of 

■ ■ ■ these second-order resonant configurations will emerge from the data, highlighting the need for 

a better understanding of the dynamical interactions between forming planetesimals. 

Subject headings: stars: fundamental parameters — stars: individual (HD41248) — stars: rotation — 
(stars:) planetary systems 

1. Introduction be in mean motion resonances (MMRs), yet even 

though such objects can be found to have period 
Resonances seem to be a common feature in na- ^.^tios that are close to known MMRs, confirma- 

ture when bodies with measurable gravitational ^j^j^ of ^j^g existence of any MMR can only be 

fields interact dynamically with one and other. j^^de by studying the system dynamically to con- 

In the Solar System, various bodies are found to ^^^ effg^^g ^^^^^ ^s hbration of the resonant an- 

gles. Examples include the 2:5 MMR of Jupiter 
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and Saturn (JMichtchenko fc Ferraz-Mellol l200l[ ) 
and the 3:2 MMR between Neptune and Pluto 
([Pealc 1976). In fact a 1:2 MMR between Jupiter 
and Saturn may have been the driving force be- 
hind the curre nt configuration of the solar systems 
outer pl anets (IGomes et al.l l2005t iTsiganis et al. 
En0,4 iMorbidelh et alJl2005h . 



Beyond planetary bodies there are also various 
moons and asteroids that exhibit MMRs. A fa- 
mous example of these comes from Jupitcrs moons 
Ganymedej_Europa, and lo that form a 1:2:4 res- 
onant set (jPeale 19 761 . Such resonances provide 
important constraints on the formation and evolu- 
tion of migrating bodies and allow a window into 
the dynamics of evolving systems. Therefore, dis- 
covering new resonant exoplanetary systems can 
provide unique constraints on the early evolution 
of exoplanets, along with a contextual view of the 
whole ensemble of planet formation and evolution 
in general. 

Currently we have confirmed a number of exo- 
planetary systems that have pairs or more of plan- 
ets in some resonant configuration by radial veloc- 
ities. 2:1 resonances seem to be the most com- 
mon by-product of giant planet formation (e.g. 



Marcv et al.ll2001l ). Further resonances have been 



witnessed, including a Jupiter-moon lik e 1:2:4 or- 



bital resonant set ( Rivera et al.l |2010|), however 



these systems tend to be first-order resonances. 

(|2008l) reported the discovery 



Desort et al 



of the HD60532 exoplanetary system that con- 
tains a pair of planets in a second-order MMR. 
Again these planets were found to be gas gi- 
ants with masse s abo ve a Jupiter-mass and 
Laskar fc Correial ( 20091 ) confirmed they are in- 
deed in a 3:1 MMR configuration. Recently, 



Fabrvckv et al.l (|2012f) used Kepler transit timing 



measurements to detect resonant systems includ- 
ing that of Kepler- 29 that seems to have two plan- 
ets orbiting the star locked in a 9:7 MMR. This 
system does not have radial velocity information 
to confirm the nature of the objects, however this 
could be the first super-Earth planetary system in 
such a dynamical commensurability. Kepler data 
also indicate that 3:1 MM Rs are the most corn mon 
second-order resonances (JLissauer et al.ll201lh . 

In this work we present the discovery of a pos- 
sible new second-order MMR configuration that 
has never been previously observed for planets. 
We show that the magnetically inactive G2V star 



HD41248 hosts at least two rocky planets that are 
on the edge of a 7:5 MMR, by analysing exten- 
sive archival data from HARPS that previously 
did not contain any known Doppler signals. In §[5] 
we present the data we have used in this study, 
then in §[3] we discuss the comparison between this 
star and the Sun. In §[4] we present the Keplerian 
analysis of our radial velocities, and in §[5] we show 
that the signals are not associated with any activ- 
ity related phenomena. Finally, § [5] presents our 
dynamical stability tests for this system and we 
finish with a discussion of the system in § [T] 

2. HD41248 Data 

All data in this paper were taken from the ESO 
HARPS Archiveo, a community tool that allows 
users to download fully reduced and analysed data 
that has been processed using the HARPS-DRS 
Version 3.5. The pipeline performs the usual re- 
duction steps for high resolution echelle spectra, 
from bias and flatfielding, to extraction and wave- 
length calibration. 

A total of 62 velocities for the star HD41248 
were downloaded and used in this analysis as part 
of our project to discover new rocky planets orbit- 
ing nearby Sun-like stars using our novel method- 
ologies. The baseline of observations is close to 
7.5 years (BJD 2452943.85284 to 2455647.57967) 
and in general a high level of data quality was 
maintained. The median S/N for the set is ~100 
at a wavelength of around 6050A, with a lowest 
value of 26 and a highest value of 150. None of 
the data we downloaded were rejected from our 
final analysis. 

After reduction and extraction of the observa- 
tion has finished, post-reduction analysis is also 
performed on the spectra and this consists of cross- 
correlating each of the echelle orders with a pre- 
fabricated binary mask to generate an order-per- 
order cross-correlation function (CCF), which are 
then combined using a weighting scheme to pro- 
duce a single stable mean CCF. This mean CCF 
is then fit by a gaussian function and the gaussian 
model allows the software to generate a precise 
and absolute radial velocity measurement. 
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Table 1 
harps radial velocities for hd41248. 



BJD 


DRS-RV 


DRS-cTRv 


TERRA-RV 


TERRA-CTRv 


Smw 


logi^HK 


BIS 


FWHM 


[days] 


[ms-i] 


[ms-i] 


[ms-i] 


[ms-i] 


[dcx] 


[dex] 


[ms-i] 


[ms-i] 


2452943.8528426 


3526.591 


2.588 


-3.512 


2.591 


0.1686±0.0024 


-4.9348 


35.92±3.66 


6721.78 


2452989.7102293 


3519.139 


4.063 


-7.736 


4.223 


0.1566±0.0045 


-5.0089 


27.39±5.74 


6719.01 


2452998.6898180 


3526.427 


5.425 


1.179 


6.226 


0.1495±0.0062 


-5.0599 


33.53±7.67 


6701.20 


2453007.6786518 


3526.626 


2.525 


-3.456 


2.232 


0.1609±0.0025 


-4.9810 


28.61±3.57 


6718.19 


2453787.6079555 


3522.442 


2.758 


-3.192 


3.353 


0.1566±0.0031 


-5.0091 


31.30±3.90 


6718.53 


2454055.8375443 


3523.175 


2.062 


-6.636 


2.226 


0.1674±0.0018 


-4.9418 


23.94±2.91 


6714.51 


2454789.7207967 


3522.987 


0.817 


-4.361 


0.897 


0.1705±0.0005 


-4.9241 


27.43±1.15 


6722.19 


2454790.6943362 


3519.487 


0.899 


-6.817 


0.916 


0.1698±0.0005 


-4.9281 


30.82±1.27 


6724.20 


2454791.7055725 


3522.466 


0.834 


-4.765 


0.831 


0.1708±0.0005 


-4.9225 


29.54±1.18 


6720.59 


2454792.7042506 


3522.290 


0.795 


-4.763 


0.779 


0.1723±0.0005 


-4.9143 


28.09±1.12 


6728.64 


2454793.7211230 


3524.992 


0.890 


-2.854 


0.984 


0.1727±0.0005 


-4.9124 


25.27±1.25 


6727.73 


2454794.6946036 


3527.038 


0.893 


0.860 


1.051 


0.1715±0.0005 


-4.9191 


29.58±1.26 


6732.04 


2454795.7156306 


3528.445 


0.913 


0.874 


0.909 


0.1737±0.0005 


-4.9070 


29.54±1.29 


6725.91 


2454796.7195391 


3528.207 


0.957 


0.513 


0.975 


0.1736±0.0006 


-4.9076 


30.32±1.35 


6727.87 


2454797.7051254 


3528.994 


0.908 


3.036 


1.137 


0.1750±0.0005 


-4.9001 


27.76±1.28 


6733.21 


2454798.6972277 


3531.198 


0.915 


3.831 


0.906 


0.1731±0.0005 


-4.9101 


25.13±1.29 


6731.22 


2454902.5907553 


3525.084 


2.021 


0.000 


2.525 


0.1796±0.0018 


-4.8768 


21.34±2.85 


6729.22 


2454903.5172666 


3527.587 


0.783 


0.359 


0.909 


0.1722±0.0004 


-4.9151 


27.68±1.10 


6726.66 


2454904.5185682 


3525.760 


0.901 


-2.032 


0.862 


0.1737±0.0005 


-4.9071 


27.18±1.27 


6722.67 


2454905.5355291 


3527.552 


0.898 


-0.148 


0.897 


0.1706±0.0005 


-4.9240 


27.97±1.27 


6722.94 


2454906.5179999 


3527.947 


1.041 


0.122 


1.114 


0.1713±0.0006 


-4.9201 


31.02±1.47 


6723.85 


2454907.5647983 


3527.517 


0.905 


0.140 


1.123 


0.1731±0.0006 


-4.9101 


28.88±1.28 


6727.04 


2454908.5603822 


3526.226 


0.804 


-2.046 


0.961 


0.1721±0.0004 


-4.9155 


27.94±1.13 


6720.84 


2454909.5380036 


3527.107 


0.854 


-0.513 


0.928 


0.1699±0.0005 


-4.9277 


25.45±1.20 


6721.35 


2454910.5385064 


3528.236 


1.115 


1.119 


1.170 


0.1727±0.0008 


-4.9121 


27.50±1.57 


6726.82 


2454911.5427244 


3524.975 


0.747 


-3.570 


0.722 


0.1715±0.0004 


-4.9189 


28.13±1.05 


6721.04 


2454912.5392921 


3525.216 


0.739 


-3.038 


0.898 


0.1719±0.0004 


-4.9167 


25.78±1.04 


6718.78 


2455284.5272133 


3528.228 


0.734 


-1.092 


0.917 


0.1751±0.0004 


-4.8998 


28.42±1.03 


6724.46 


2455287.5109103 


3524.603 


0.881 


-4.445 


0.921 


0.1743±0.0005 


-4.9038 


25.18±1.24 


6733.01 


2455288.5285775 


3523.316 


0.754 


-4.407 


0.807 


0.1731±0.0004 


-4.9104 


25.43±1.06 


6730.40 


2455289.5460248 


3526.698 


0.789 


-1.501 


0.923 


0.1743±0.0004 


-4.9036 


26.29±1.11 


6723.59 


2455290.5095380 


3525.900 


0.891 


-0.647 


0.996 


0.1737±0.0005 


-4.9068 


21.95±1.26 


6727.32 


2455291.5216615 


3526.083 


1.006 


-3.054 


1.060 


0.1707±0.0006 


-4.9233 


26.85±1.42 


6734.67 


2455293.5043818 


3527.359 


0.974 


1.254 


1.113 


0.1719±0.0006 


-4.9165 


27.39±1.37 


6727.22 


2455304.5180173 


3522.412 


1.689 


-0.532 


5.028 


0.1318±0.0038 


-5.2203 


31.20±2.38 


6801.22 


2455328.4550220 


3529.109 


0.791 


2.304 


0.828 


0.1746±0.0005 


-4.9022 


20.73 ±1.11 


6735.38 


2455334.4564390 


3532.003 


1.365 


7.425 


1.432 


0.1685±0.0014 


-4.9353 


25.33 ±1.93 


6737.32 


2455387.9305071 


3530.983 


1.039 


4.143 


0.986 


0.1708±0.0008 


-4.9225 


27.29 ±1.47 


6737.65 


2455390.9312188 


3531.330 


1.577 


4.830 


1.517 


0.1635±0.0013 


-4.9652 


30.24 ±2.23 


6734.97 


2455434.8790630 


3516.749 


3.334 


-4.089 


3.355 


0.1479±0.0036 


-5.0724 


31.59±4.71 


6740.53 


2455439.8843407 


3527.755 


3.258 


10.200 


4.319 


0.1488±0.0038 


-5.0649 


32.80±4.60 


6723.47 


2455445.9241644 


3522.932 


3.113 


-0.217 


3.839 


0.1545±0.0036 


-5.0233 


28.30±4.40 


6727.79 


2455465.8566225 


3526.615 


1.412 


1.415 


1.440 


0.1706±0.0012 


-4.9237 


23.18 ±1.99 


6731.85 


2455480.8795598 


3527.592 


1.133 


1.534 


1.245 


0.1706±0.0008 


-4.9237 


27.09 ±1.60 


6728.47 


2455483.8136024 


3525.010 


1.740 


-1.948 


1.837 


0.1704±0.0015 


-4.9250 


27.97±2.46 


6734.01 


2455488.8262312 


3525.839 


0.747 


-0.835 


0.818 


0.1708±0.0004 


-4.9226 


31.07±1.05 


6725.47 


2455494.8532009 


3529.808 


0.954 


3.114 


1.123 


0.1686±0.0006 


-4.9348 


26.88 ±1.34 


6730.94 


2455513.7823100 


3528.860 


1.286 


-0.462 


1.438 


0.1715±0.0009 


-4.9186 


35.60±1.81 


6740.83 


2455516.7515789 


3530.744 


0.936 


4.112 


1.032 


0.1749±0.0007 


-4.9005 


22.71 ±1.32 


6736.78 


2455519.7046533 


3528.854 


1.198 


3.453 


1.116 


0.1682±0.0009 


-4.9375 


21.95 ±1.69 


6735.81 


2455537.7997291 


3527.714 


0.718 


0.897 


0.859 


0.1756±0.0005 


-4.8968 


27.01 ±1.01 


6730.29 


2455545.7213656 


3529.198 


0.854 


3.099 


0.990 


0.1749±0.0005 


-4.9006 


26.68 ±1.20 


6737.59 


2455549.7548559 


3527.949 


0.833 


-0.567 


0.896 


0.1714±0.0005 


-4.9195 


31.97±1.17 


6734.54 


2455576.7923481 


3525.346 


1.014 


-1.481 


1.120 


0.1654±0.0007 


-4.9533 


29.17±1.43 


6733.68 


2455580.7312518 


3519.411 


0.895 


-7.373 


0.965 


0.1688±0.0006 


-4.9339 


32.94±1.26 


6729.42 


2455589.7734088 


3528.423 


1.420 


0.958 


1.579 


0.1684±0.0013 


-4.9360 


24.59 ±2.00 


6733.73 


2455612.6068850 


3527.663 


0.816 


0.976 


0.713 


0.1714±0.0005 


-4.9193 


25.58 ±1.15 


6731.09 



The values extracted from the ESO Archive 
DRS processing are also backed up with the veloci- 
ties generated using the HAR PS-TER RA software 
(lAnglada-Escude fc Butleij2012HAnglada-Escude et al 



2012I ). which we use as a sanity check for stars like 
the Sun as the DRS values tend to produce higher 
precision for these stars, in comparison to cooler 
M dwarfs where TERRA works better. The fi- 
nal Barycentric Julian Dates, DRS, and TERRA 
radial velocities, along with DRS and TERRA 
uncertainties, are shown in Table [2l 

3. HD41248 vs the Sun 

The properties of HD41248 are summarised in 
Table [5] but some of the most interesting features 
are highlighted here. First of all, the star has a 
spectral type of G2V and B — V colour of 0.624, 
meaning it is a solar analogue since the solar val- 
ues are G2V and 0.642, respectively. Rocky plan- 
ets around such stars can provide direct tests of 
planet formation mechanisms and architectures of 
systems orbiting Sun-like stars. 

3.1. Chromospheric Activity 

Low amplitude signals are most easily sought 
after in radial velocity data of the most quies- 
cent and slowly rotating stars. HD41248 has a 
logi?HK activity index of -4.92 dex and a rota- 
tional velocity of only 2.4± kms~^, highlighting 
that this star is an ideal candidate for such studies. 
The Sun has a logi?HK acti vity index of -4.91 dex 
and vsini of 1.6±0.3 kms~^ (|Pavlenko et al.ll2012n 
showing both these stars have good agreement in 
their evolutionary properties too. 

Fig. [T] shows the distribution of these activ- 
ity indices and the best fit model Gaussian to 
the data. The log-RuK values were computed us- 
ing the HARPS ID spectra follow i ng the pro- 
cedur es explained in Ijenkins et al.l (J2006l l2008l 
20111 ). Clearly most of the data are tightly clus- 
tered around the mean logi?HK of -4.92 dex. The 
distribution is tightly packed, with a scatter of 
only 0.01 dex, which is both intrinsic variability 
of the star and uncertainty in the measurements. 
There are a few values with activities less than - 
5 dex, however these tend to be the lowest S/N 
spectra and therefore there is a bias in the line 
core measurements as random noise causes nega- 
tive read values that artificially draw the line core 



values lower than they really are, and therefore the 
overall activity index is found to be lower than it 
should be. 

3.2. Abundance Pattern 

One parameter where there is some difference 
between the solar value and that of HD41248 is the 
stellar metallicity. The iron abundance ([Fe/H]) of 
this star is -0.43±0.10 dex, which is significantly 
lower than the solar value, defined as 0.0 dex. In 
fact, this low abundance value could be one of 
the important parameters that defines a low-mass 
rocky dominated system from a gas giant dom- 
inated system, since the rising power law prob- 
ability function that shows the most metal-rich 
Sun- like stars have a higher probability of hosting 
giant planets (JFischer fc Valentill2005l : ISousa et al 



20111 ) seems to disappear o r turnover for the rocky 



super-Earth popul ation (jBuchhave et al.l 12012 : 



Jenkins et al 



20131 ). Therefore, HD41248 can be 



classed as a metal-deficient and old solar analogue. 

Table [3] lists various abundances for volatile el- 
ements in the atmosphere of HD41248. We mea- 
sured these values directly from the spectra us- 
ing the method explained in lPavlenko et al.l ( 2012 ) 
and more details of how we arrived at the se values 
will be discussed in llvanyuk et al.l (20131). These 
abundances generally track the low [Fe/H] values, 
and in comparison to the Sun, all of the elements 
we have considered here are depleted. This shows 
that the nascent disk from which planets could 
form was metal-deficient, are therefore depleted in 
the typical elements we expect are processed into 
planetesimals to form systems of planets through 
core accretion. 

Finally, it is necessary to know the stellar 
mass of HD41248 so that we have a handle on 
the masses of any planets detected orbiting the 
star from the Doppler curve. Given the Te// we 
measure of 5713±50 K, the absolute magni- 
tude (Mv) of 5.22±0.08 mags, and the metal- 
licity of -0.43±0.1 dex, we find a stellar mass 
of 0.92±0.05 Mq. This is in agreement with 
other work s who find s i milar values, both lower 
(0.81 M(7); [S ousa et all l201l[ ) and higher (0.97; 



Casagrande et al. 20111 ). 



Table 1 — Continued 



BJD 

[days] 



DRS-RV 



D RS -cTftv 



TERRA-RV 

fms-M 



TERRA-cTRv 

fms-i] 



Smw 
[dox] 



logi?,HK 

[dox] 



BIS 



FWHM 

fms-M 



2455623.6361828 


3528.301 


1.175 


0.455 


1.316 


0.1685±0.0011 


-4.9357 


28.96 ±1.66 


6726.23 


2455629.5528393 


3528.652 


0.944 


1.864 


1.002 


0.1682±0.0006 


-4.9370 


25.71 ±1.33 


6732.34 


2455641.5542106 


3528.979 


1.005 


1.879 


1.012 


0.1697±0.0007 


-4.9286 


28.91 ±1.42 


6727.20 


2455644.5845033 


3529.351 


1.260 


3.245 


1.381 


0.1742±0.0006 


-4.9043 


19.09 ±1.78 


6740.10 


2455647.5796694 


3529.945 


1.372 


1.898 


1.259 


0.1669±0.0012 


-4.9448 


29.10 ±1.94 


6737.56 
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Fig. 1. — Histogram of the logi?HK activity indices 
of the timeseries data for HD41248. The best fit 
Gaussian model to the activities is represented by 
the sohd curve and the associated data for this 
model is shown in the plot. 



Table 2: Stellar parameters for HD41248. 



Parameter 


HD41248 


Reference 


RA J2000 (h:m:s) 


06" 00"" 32=. 781 


Perrvman et al. C 19971 


Dec J2000 (d:m:s) 


-56°09'42.61" 


Perrvman et al. (19971 


Spectral Type 


G2V 


Perrvman et al. (19971 


B -V 


0.624 


Perrvman et al. ^^971 


V 


8.82 


Perrvman et al. ('l9971 


distance (pc) 


52.38±1.95 


van Lccuwcn (20071 


Mv 


5.22±0.08 


This Work 


logR'uK 


-4.94 


This Work 


Hipparcos A^obs 


106 


Perrvman et al. (19971 


Hipparcos a- 


0.015 


Perrvman et al. ^^971 


AMv 


-0.700 
0.68±0.03 


Jenkins et al. (2011^ 


LJLq 


This Work 


MJMq 


0.92±0.05 


This Work 


RJRe 


0.78±0.04 


This Work 


Tetf (K) 


5713±50 


vanvuk et al. (20131 


[Fe/H] 


-0.43±0.10 


vanvuk et al. (20131 


log g 


4.48±0.10 


vanvuk et al. (20131 


U,V,W (km/s) 


-12.97, -17.81, 26.74 
16 


Jenkins et al. (2011) 


frot.v Bin i (days) 




Dsin2(km/s) 


2.4±0.2 


Ivanvuk et al. (20131 


Age (Gyrs) 


2tl 


This Work 


Jitter - fit (m/s) 


0-90tlit 


This Work 



4. Doppler Analysis 

The radial velocity timeseries of 62 observa- 
tions for HD41248 show evidence for at least 
two low-amplitude signals embedded in the data. 
Both a periodogram analysis and our Bayesian fit- 
ting method detected these signals, however the 
Bayesian method can detect the second signal with 
a high degree of significance such that we can con- 
firm the signal is robust. 

4.1. Periodogram Analysis 

A periodogram search for strong and stable fre- 
quencies in the radial velocity dataset for HD41248 
reveals a significant peak around 18 days. The 
top panel of Fig. [2] shows these frequencies, with 
the horizontal dashed line marking the 1% false- 
alarm probability (FAP) and the horizontal dot- 
dashed line marking the 0.1% FAP, both were 
mea sured using a bootstrap an alysis technique 
fsee lAnglada-Escude et al.l 120121 ). It can be seen 
that the signal is stronger than both these bound- 
aries. 

The FAPs from the bootstrap analysis allow 
us to understand the significance of the frequency 
peaks we detect in the data. Bootstrapping is ap- 
propriate in this case since we can generate test 
data sets from the original data set, without as- 
suming any underlying distribution for the veloc- 
ities, or more importantly, their uncertainties. A 
lot of effort is being spent at the present time to 
understand how the combination of white and red 
noise affects the overall uncertainty we can assume 
for any radial velocity time series, particularly in 
HARPS data (JBaluevI IJOll tuomilSb), however 
it is still an extremely difficult task to model such 
data with any high degree of accuracy. 

We resampled the HD41248 data lO'OOO times 
with replacement, using a Monte Carlo approach 
to scramble the timestamps of the velocities, but 



retaining the velocities and their associated uncer- 
tainties. With each new sample, we recompute the 
LS periodogram and measure the strength of the 
strongest peak that we find. The strength of this 
strongest peak is then compared to the strength of 
the original peak from the observed data set and 
the FAP relates to the number of times a peak 
stronger than this one is found. In this way we can 
directly measure the probability from the data in 
an unparametric way. 

The center panel in Fig. [2] shows the peri- 
odogram of the residuals to the best fit Keple- 
rian of ~18 days. There is a clear emerging sig- 
nal around 25 days, indicating another physical 
process is causing a frequency peak at this pe- 
riod, however the significance is below the 1% FAP 
level, meaning the archival data we have used is 
not yet abundant enough to confirm the nature of 
this peak. 

The lower plot in Fig. [3] shows the periodogram 
power for the widow function, with the strongest 
energies found at long periods, beyond the baseline 
of the data. No sampling power is found around 
18 days or 25 days, indicating the periods we de- 
tect in the radial velocities are not sampling fea- 
tures. Also, no sampling power was found around 
65 days, which is an alias that could arise since 
g^ « yJ-^ — 2^ J however we do find some ac- 
tivity power around 60-70 days that we discuss 
later but our Bayesian analysis indicates this is 
not the source of the secondary signal and it is a 
real Doppler shift, not an alias. Such a tantilising 
system necessitates another methodology to con- 



Table 3: Chemical abundances for HD41248. 



Element 


[X/H] 


[X/H]0 


[X/Fe] 




[dex] 


[dex] 


[dex] 


Si I 


-4.743±0.031 


-0.34 


0.09 


Si 11 


-4.595±0.025 


-0.20 


0.24^ 


Cal 


-5.972±0.017 


-0.31 


-1.14 


Til 


-7.419±0.017 


-0.44 


-2.59^ 


Tin 


-7.326±0.032 


-0.35 


-2.49 


Vi 


-8.348±0.216 


-0.31 


-3.52^ 


Cri 


-6.791±0.020 


-0.35 


-1.96 


Pel 


-4.833±0.007 


-0.42 


- 


Fe 11 


-4.845±0.018 


-0.44 


-0.02 


Nil 


-6.206±0.010 


-0.39 


-1.37 
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Fig. 2. — Lomb-Scargle periodograms of the 
HD41248 radial velocity dataset is shown in the 
top plot, with the same analysis performed on the 
residuals to the 1-planet fit shown in the lower 
plot. The dashed and dot-dashed lines mark the 
1% and 0.1% FAP limits. The lower plot shows 
the window function, along with two vertical dash 
lines that show the positions of the two signals. 



firm the significance of this signal and we turn to 
our Bayesian analysis method since this has been 
shown to efficiently detect significant signals with 
less velocity rich data than can be accomplished 
usin g a standard Lomb-Scar gle periodogram anal- 
ysis (ITuomi fc Jenkind |2012|) . 



4.2. Bayesian Search 

Bayesian signal detection techniques have re- 
cently been introduced to spearhead the detection 
of radial velocity signatures of low-mass exoplan- 
ets around nearby stars (e.g . lAngl ada-Escude fc Tuomil 
2OI2I: iTuomi fc Jo"^ 12012 : Tuo mi et al.ll2013bf ). 
yet they have been shown to be rather immu ne to 
detections of false positives (e.g. lTuomill201 II ) . We 
analysed the HD41248 velocities using posterior 
sampling techniques and Bayesian model selection 
to find the best statistical descriptions, i.e. mod- 
els, of the data, and to obtain estimates for the pa- 
rameter probabilit y densities of the corre sponding 
models. Following iTuomi fc JonesI (|2012r ). we per- 
formed the sa mplings using the ad aptive Metropo- 
lis algorithm (JHaario et al.l 120011 ) and calculated 
the Bayesian evidences of each model using their 
truncated posterior mixture estimates. As we 
analysed the data in the Bayesian context, we 
defined the prior probability densities and model 
proba biliti es according t o the choices of iTuomi 
(|2012l) and lTuomierall ()2013bl) . 

When modeling the data, we used a common 
Gaussian white noise model as a reference model 
and attempted to improve this description by in- 
cluding correlations within it. This means that the 
measurement mean is described using a function 
IJ,{9, ^, i), where 9 is the parameter vector, t repre- 
sents time, and ^ is a vector containing any other 
variables that might have an effect on the func- 
tion. With this notation, we defined this model 

as 

p 



M(0,c,o = /fe(^p,o+7+ Vc,e 



(1) 



where fk {t) is a function describing the superposi- 
tion of k Keplerian curves with orbital parameters 
9p, 7 is a reference velocity w.r.t. the data mean, 
and parameters Ci describe the linear dependence 
of the function on the variables ^i. We used the 
S-index, BIS, and FWHM as these variables to be 
able to take into account the correlations of the ra- 
dial velocities with these activity-related indices. 



At this point we can introduce the activity indi- 
cators that we use in our Bayesian model. The BIS 
values and the FWHMs are taken directly from 
the HARPS-DRS output and details of their origin 
and usefulness can be fou nd in lOueloz et al.l (|2001l ) 



and Santos et al. 



(|2010[) . The chromospheric 5- 
indices have been measured following our own 
recipes, as mentioned in § [3l Using this proce- 
dure the uncertainties of these HARPS S'-indices 
are found to be less than 1%. 

Our results indicate that there are two signifi- 
cant periodicities in the HD41248 velocities at 18.4 
and 25.6 days. We demonstrate the significance of 
the two signals by showing the log-Bayesian ev- 
idences of models with up to two Keplerian sig- 
nals (Table |4]). Accounting for the correlations 
between the velocities and the activity indices im- 
proves the model clearly for fc = and fc = 1 but 
only marginally for fc = 2. Yet, the two-Keplerian 
model is clearly the best description of the data 
regardless of whether we account for these corre- 
lations or not and the second signal is detec t ed ac- 
cording to the detection criteria of iTuomil (J2012l ) 
because the two-Keplerian model is 3.3 x 10"^ times 
more probable than the one-Keplerian model. The 
other two detection criteria, i.e. that the orbital 
periods and radial velocity amplitudes are well- 
constrained, are also satisfied. We demonstrate 
this by plotting the estimated posterior densities 
of the orbital periods, velocity amplitudes, and or- 
bital eccentricities in Fig. [31 We have also tabu- 
lated our estimates for model parameters in Ta- 
ble H 

Obtaining the samples from the posterior den- 
sities was simple in this case because the period- 
space contained only two significant maxima cor- 



Table 4: Model probabilities 

~k hTBE hTBE P(fc) P(fc) 

(Ref) (Cor) (Ref) (Cor) 

~0 -158.1 -154.9 5.3 xlO"^^ 6.0 xlO"^^ 

1 -141.5 -136.5 4.5 xlO^^ 3.0 xlO"* 

2 -128.5 -127.7 - 1 - 1 



Note. — Log-Bayesian evidences and associated probabilities 
for the Keplerian models from k = 0, ..., 2 for the HD41248 
velocities taking into account the correlation terms with activity 
indices (Cor) and without including these terms (Ref). 
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Fig. 3. — The left column corresponds to the short period signal and the right column corresponds to the 
longer period signal. The top two histograms show the estimated probability density function of periods 
from the Bayesian analysis, the middle histograms are the same for the semi-amplitude of the signals, and 
the lower two histograms are for the eccentricities of the signals. Also overplotted by the solid curves are 
Gaussians with the same mean and sigma as the distributions. 



Table 5: Keplerian solutions for HD41248. 



Paramctcr 


HD41248 h 


HD41248 c 


P [days] 


18.357 [18.313, 18.423] 


25.648 [25.518, 25.768] 


e 


0.15 [0, 0.26] 


0.00 [0, 0.18] 


w [rad] 


3.2 [0, 27r] 


1.0 [0, 27r]i 


Afo [rad] 


0.2 [0, 27r] 


1.7 [0, 27r]i 


K [ms-i] 


2.93 [1.65, 3.65] 


1.84 [0.67, 2.97] 


a [AU] 


0.137 [0.126, 0.154] 


0.172 [0.158, 0.192] 


rUp sinz [M0] 


12.3 [6.9, 16.5] 


8.6 [3.6, 15.1] 


7 [ms"-^] 


-48 [-1.36, 0.18] 




CTj [ms^i] 


0.90 [0.57, 1.84] 




ci [ms^-^ dex^-^] 


147 [-20, 296] 




C2 


0.000 [-0.232, 0.229] 




C3 [10-3] 


61 [-14, 156] 






5 10 15 

Orbital phase, signal 01 [days] 



Note. — The maximum a posteriori estimates of the solution for HD41248 velocities and the 99% Bayesian credibility intervals 
from posterior samplings. 



responding to the two signals we observed. De- 
spite several attempts, we could not find a third 
signal and the samplings of the parameter space 
of a three-Keplerian model did not converge to a 
third periodicity. We also attempted to include a 
moving average (MA) component in the statisti- 
cal model to improve its performance and to take 
into account intrinsic correlation s in the measure- 
ment noise (JTuomi et al.ll2013bl lal). However, the 
corresponding parameter describing the amount of 
autocorrelation in the data was found to be con- 
sistent with negligible estimates, and the corre- 
sponding MA component did not improve the sta- 
tistical model. Furthermore, as can be seen in 
Table O all the parameters quantifying the linear 
dependence of the velocities on the different activ- 
ity indices (ci, C2, and C3) are consistent with zero 
(i.e. cannot be shown different from zero with a 
99% credibility), which explains why taking these 
correlations into account improved the statistical 
model little in terms of Bayesian evidences (Ta- 
ble 2]). However, there are likely correlations be- 
tween the velocities and S-index and FWHM, but 
BIS cannot be shown to be correlated with them 
at all. 

Based on the Bayesian analyses of the HD41248 
radial velocities, there are two significant periodic- 
ities in the data corresponding to two super-Earths 
or Neptune-mass planets on nearby close-in orbits 
(Fig. 2]). In fact, their orbits arc so close to one 
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Fig. 4. — Phased radial velocity curves for plan- 
ets HD41248 h and c from the Bayesian detected 
periods. 



another that the 99% credibility intervals of Of, 
and Oc taking into account the uncertainty in the 
stellar mass are almost overlapping (Table [5|) and 
the configuration cannot be immediately stated to 
be stable in the long term. However, the orbital 
stability is suggested by the fact that the ratio 
of the orbital periods is almost exactly 1.4 with 
a Bayesian 99% interval of [1.384, 1.407], which 
suggests a possible 7:5 MMR. The eccentricities 
of the system also show an interesting configura- 
tion with the inner planet having an eccentricity 
around 0.2 and the outer planet having an eccen- 
tricity close to circular. In fact, we tested if the 
inner planet's eccentricity was closer to zero by 
changing the eccentricity prior model to favour a 
more circular orbit, but the eccentricity was found 
to be 0.22 in this case, indicating the inner planet 
does have some genuine measurable eccentricity 
given the data. If we can rule out the source of 
these signals as originating from stellar magnetic 
activity, this could be a remarkable 7:5 MMR sys- 
tem, and therefore we assess this possibility in § |6] 
by analysing the stability of the system. 

5. Activity Indicators 

Within our Bayesian model we have taken into 
account any linear correlations between the ra- 
dial velocities and the activity indicators (BIS, 
FWHM, and 5harps): as explained in the previ- 
ous section. However, not only does the presence 
of activity add additional noise to any individual 
radial velocity measurement, but it can also mimic 
the presence of a true Doppler induced frequency 
in the full timeser ies of radial velocity data (e.g. 



Queloz et al.l 120011 ). Therefore, it is necessary to 



perform a test on whatever activity indicators one 
has at hand, to search for periodicities that could 
be associated with the periodicities found in the 
velocities. 

In Fig. [5] we show the periodograms for each of 
these activity indices, going from the S'-indices at 
the top, to the BIS values in the center, and finally 
the FWHMs at the bottom. The 5-indices show a 
strong periodogram peak at 61 days and another 
slightly weaker peak at 69 days which could be 
the rotational period of the star, or some other 
activity cycle. There is no significantly strong fre- 
quency peak in the BIS periodogram, however in 
the FWHM periodogram there is a strong peak 
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Fig. 5. — Periodograms of the activity indices we 
consider. From top to bottom we show the S'- 
indices, BIS values, and the FWHM values. 
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emerging at ^-^71 days that could be related to the 
same frequency detected in the S'-indices. If so 
this may back-up the hypothesis that the rota- 
tional period of this star, or a harmonic thereof, 
is around 70 days, or at least there is an activity 
cycle present at that frequency. 

The inner signal in the radial velocities at 
18 days appears close to a peak in the BIS peri- 
odogram at around 18 days and this could throw 
some doubt on the nature of the inner signal, even 
though it is very unlikely for two signals to appear 
in a radial velocity dataset entrenched so close to 
such a MMR, down to the level of minutes. Also 
the peak in the BIS periodogram is not significant, 
with many more stronger peaks in the data, and 
since we included correlations between the BIS 
and the radial velocities, we can be confident this 
is not the source of the short period signal we de- 
tect. We also note that any power in the activity 
indicators at a frequency that could give rise to an 
alias in a periodogram search for signals, will not 
be detected by our Bayesian analysis, since our se- 
lection must give rise to signals with amplitudes 
that are significantly different from zero. This also 
protects our method from being confused by peaks 
that we see in the window function. Therefore, 
we must conclude that the two signals we detect 
close to a 7:5 MMR are genuine Doppler signals, 
induced by the presence of two low-mass planets 
orbiting HD41248. 

6. System Stability 

The HD41248 system poses an interesting dy- 
namical challenge because of the proximity of the 
7:5 mean-motion resonance and the close spacing 
of the two planets, both to each other, and to the 
central star. Here General Relativistic (GR) ef- 
fects play a role in the dynamical stability. Before 
moving on to perform a numerical search for sta- 
ble configurations, we list several key properties of 
the system that determine its stability. 

Chambers! (|200ll) defined several measures that 



he used to quantify a system of planets, th ese are: 
the Angular Momentum Deficit (AMD; iLaskar 
19971 ). the fraction of total mass in the most mas- 
sive planet (Sm), a spacing parameter (5s) that 
scales as the planet to star mass ratio fj,^^'^ rather 
than the Hill relation of fj}^^ ( Chambers et al. 
1996I ). and a concentration parameter {Sc) which 



measures how the mass is concentrated in an annu- 
lus. We supplement this with the measure stating 
the average spacing in Hill radii (Sh)- The values 
of these quantities for Venus and Earth, Jupiter 
and Saturn and HD41248 are hsted in Table i] 

As one may see, the AMD of HD41248 is nearly 
twice as high as that of Jupiter and Saturn, and 
much higher than that of Venus and Earth. Sys- 
tems with hig her AMD ha ve the possibility to be 
more chaotic (|Laskaiill997l ) and have more oppor- 
tunity to exchange it among the planets. The rel- 
atively high AMD of HD41248 and the fact that 
one planet appears much more eccentric than the 
other suggests that both eccentricities will fluctu- 
ate with a large amplitude, with one planet being 
at a minimum (c) when the other is at a maximum 
(b). 

Given that both planets have a nearly equal 
mass it is no surprise that Sd ^ 0.5. More inter- 
esting is the spacing parameter. This quantity is 
only ^7 while for Venus and Earth it is about 18. 
Even for Jupiter and Saturn it is slightly larger. 
However, it is interesting to see how Sh scales as 
a function of the masses of the planets. We have 
Sh = {dc — o,b)/fH where th is the mutual Hill 
radius i.e. r^ — ^(a^ + ac)[{mb + mc)/3A/*]"'^/'^. 
This can be solved for ab/ac to give 



ttr 



2T-S 



H 



2T + Sh 



(2) 



where T = [3M*/(mh-|-mc)]^/^. Setting Sh = 7 
and TTifc ^ rric '^ 10 M^ yields ab/ac ^ 0.83, but 
when inserting the masses of Jupiter and Saturn 
one has aj /as ~ 0.58. Thus, even though Jupiter 
and Saturn have a similar spacing than HD41248 
in terms of their mutual Hill radii, they are spaced 
father apart in relative semi-major axis ratio than 
the HD41248 system, implying that Jupiter and 
Saturn need a higher eccentricity to begin cross- 
ing their orbits than the planets of HD41248. In- 



System 


AMD 


^7n 


Ss 


Sc 


Sh 


VE 

JS 

HD41248 


8 xlO"4 
1.6 xlO-^ 
3.4 xlO-3 


0.554 
0.769 
0.560 


17.5 
8.3 
6.9 


204.5 
80.3 
411.0 


26.2 
7.9 

8.5 



Table 6: IChambersI ()2001l) stability quantities for 
the terrestrial planets, giant planets and HD41248. 
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deed, for HD41248 the orbits begin to cross when 
both planets have an eccentricity near 0.1, so un- 
less they are apsidally aligned or protected by a 
resonance, they will become unstable. The nomi- 
nal eccentricity of planet b is near 0.1, so the sys- 
tem is close to instability and the stability may 
crucially depend on their phasing. 

In Fig. [6] we have plotted the secular evolu- 
tion of both planets using the nominal orbital el- 
ements. The value of Uc is unconstrained so we 
set An7 = ujh — ujc i.e. the difference in the argu- 
ment of periastron, to 20° . The right panel shows 
the evolution of planet b in the eccentricity-A-cu 
plane while the right panel shows the same for 
planet c. The planets are in apsidal alignment i.e. 
Aw librates around and both planets show sub- 
stantial excursions in their eccentricity amplitude. 
This alignment prevents many close encounters. 
We did not find the resonant angles to librate for 
this simulation. The question is whether this mo- 
tion persists for all initial conditions allowed by 
the observations. Numerical simulations should 
shed light on this issue. 

6.1. Numerical methods 

We proceeded to perform a large number [^ 
40 000) of short simulations of the HD41248 sys- 
tem, which were then analysed for stability and or 
chaos. 

We performed grid searches of both planets in 
the a — uj and a — e planes. The stability de- 
pendence upon the mutual inclination was not 
tested. At this stage the dependence on the masses 
was not tested either because of the large possible 
range depending on the orientation of the system 
with respect to the observer. However, an estimate 
of the stability dependence on the masses can be 
made and maximum masses can be computed as 
follows. 

For a planet on a fairly eccentric orbit, the 
cha otic boundary surround ing this planet is given 
bv (|Mustill fc Wvatdlioii) 



5a 



1.8(eM)^/^ 



(3) 



This differs from the /z^/'^ law of lWisdomI (|l980[ ) 
since it takes eccentricity into account. Taking 
both planets as 10 M0 and setting their eccentric- 
ities ^ 0.05 we have 5a/a ^ 0.12. The planets are 



spaced 0.035 AU apart and thus 5a/ a < 1 — ahjac-, 
suggesting the system could be stable. Chaotic 
instability would set in when 5al a ^^ 1 — ajy/a^^ 
which can be solved and yields e/i > 2 x 10~^. 
Taking e ~ 0.1 implies that the masses need to 
exceed 60 M®, well outside the uncertainties. 

We are dealing with two planets on planar or- 
bits, so that we have four degrees of freedom: the 
semi-major axis or period ratio, the eccentricities 
and An7. We performed a grid search for all of 
these parameters. 

The first grid search was done in the semi- 
major axis-argument of periastron plane of planet 
b, keeping all other orbital elements at their nomi- 
nal values. The inclination and longitude of the as- 
cending node were set to zero. For planet c Fig. [3] 
shows that Wc is almost arbitrary, so we set it to 
for the sake of simplicity. We varied Wf, evenly 
in steps of 3° between and 360° and sampled a^ 
assuming a Gaussian distribution with mean and 
standard deviation given in Fig. |3] using 99 steps. 
This resulted in 11880 simulations. The sampling 
method employed here takes into account the den- 
sity distribution of the semi-major axis but does 
not consider mutual correlations between the ele- 
ments. 

The planets were integrated for 20 500 yr using 
the S WIFT MVS integrator ()Levison k Duncan 
I994I ). We included the effects of GR by adding 



the potential term V = -3( GM Jc)^a(l - e^)/r 



where c is the speed of light (|Nobih fc Willlll986[) . 
This term reproduces the general relativistic effect 
of Mercury's perihelion advancement. For all sim- 
ulations the time step was set to 2x10""' yr and 
output was every 10 yr. 

We analyse the stability of the system using fre- 
quency analysis ( Laskar 1993; ). The basic method 
is as follows: using a numerical integration of the 
orbit over a time interval of length T, we com- 
pute a refined determination of the semi-major 
axes fli, a2 obtained over two consecutive time 
intervals of length Ti = Ta == T/2. The sta- 
bility index D — TOax(|l — 02/01!) provides a 
measure of the chaotic diffusion of the trajectory. 
Low values close to zero correspond to a regular 
solution, while high val ues are synonymous with 
strong chaotic motion (|Laskaii ) . The advantage 
of the frequency analysis method is that it does 
not require long-term simulations and thus large 
regions of phase space can be tested with a reason- 
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able amount of CPU power (jCorreia et al.ll2005l) . 
A similar methodology was employed for the 
other three sets of simulations. The eccentricities 
were sampled according to the best-fit distribu- 
tions and we set Am =10°. 

6.2. Results 

In this section we present the results from the 
numerical simulations determining the stable re- 
gions. These results will be presented in a series 
of figures. 

Fig. [7] displays the stability index of the system 
as a function of the period ratio Pc/ Pt and Atu. 
For this figure we varied the value of at within 
its observed range. The axis labels indicate the 
quantity that was varied while the colour coding 
shows logD. Motion where log I? > 10^"^ is con- 
sidered chaotic and possible instability may occur 
(jCorreiaet al.ll2005l ). 

An interesting feature are the seven horizontal 
regions of strong chaos and instability that appear 
independent of period ratio but occur at regular 
intervals as Aw varies. We verified that in these 
regions the motion is unstable: the planets en- 
counter each other on short time scales. Varying 
ttc and Wc results in essentially the same outcome 
as that displayed in Fig. [T] The structure at low 
Avj and low period ratio marks the limits of the 
7:5 resonance. 

The reason for these regions of unstable motion 
have to do with the proximity of the 7:5 resonance. 
Without the resonance the system would only have 
a single degree of freedom, Aw, and th us be in- 
tegrable (JMichtchenko fc Malhotrall2004l) . The in- 
tegrable system may exhibit one of three types 
of motion: Azu librates around (apsidal align- 
ment), Atu circulates and Aw librates around 
180° ( apsidal ant i- alignment; see iMorbidelli et al 



(|2009l ) for a discussion on these three types of mo- 
tion). Without the presence of the resonance in- 
creasing the initial value of Aw from to 360° 
results in the system displaying all three types of 
motion in three distinct regions. 

It should be noted that regions surrounding 
the period ratio of these planets are clustered 
with potential MMRs. Nearby and stronger first- 
order MMRs are the 3:2 and 4:3 commensurabili- 
ties, the similar strength 6:4 and 8:6 second-order 
MMRs around found around this dynamical re- 
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Fig. 6. — Secular evolution of the two planets of 
HD41248. The left panel shows the evolution of 
planet b in the e — Aw plane while the right panel 
shows the same for planet c. The planets are in 
apsidal alignment (Azu librates around 0). 




Fig. 7.— Stability map of the HD41248 system 
as a function of the period ratio (horizontal axis) 
and Aw (vertical axis). The scale on the right 
indicates the value of the stability index D. Blue 
regions are stable, red is unstable. We varied at 
and ujb. 



13 



gion, and so are the weaker 10:7 and 11:8 third- 
order MMRs. Therefore, we ran some tests to 
ensure the 7:5 MMR is the most hkely dynami- 
cal configuration for this system, given the orbital 
parameters. Taking the periods and uncertainties 
from Table m and assuming Gaussian distributions 
for each, we computed the distribution of the pe- 
riod ratio (Pc/Pb) and found a mean and a of 1.397 
and 0.003, respectively. The narrow range of pe- 
riods given by the orbital elements constrains the 
possible MMR to be only the 7:5 ratio, since even 
the nearest 10:7 and 11:8 MMRs are found to be 
8cr and 6a away, respectively. 

From the analysis of the motion of Atjj we 
find that the system most likely has Atn librat- 
ing around because this configuration yields all 
stable solutions. When Am circulates or librates 
around 180° the system is chaotic, and often un- 
stable. 

However, the presence of the resonance changes 
the system's behaviour. The 7:5 is second order, 
and thus there are three resonant arguments: 



(4) 



where A = nj -|- M is the mean longitude and 
vu = fl + uj and M is the mean anomaly. The 
system is planar and thus fl is undefined (we set 
it to zero) . The proximity of the resonance causes 
the seven maxima and minima in Z? as a function 
Am as it circulates or librates around 180°. 

We have decided to display this evolution in 
a series of figures. Starting from the bottom-right 
corner of Fig.[7]the system becomes unstable when 
Am ~ 30° so a transition must occur. We dis- 
played the evolution in the Cb — Am (left column) 
and Cc — Atu planes (right column) for three sim- 
ulations in Fig. m The top panels depict the evo- 
lution while the motion is still regular, the middle 
panels show the motion at the edge of the instabil- 
ity strip, and the bottom panels show the motion 
for an unstable configuration (although we only 
depicted the motion until the instability occurred). 

In the top panels both planets have Am librate 
around 0. Increasing the initial value of Am from 
the top panels to the middle panels one sees a fun- 
damental change in the middle-left panel. There 
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appear to be two regions of motion: one in which 
Am librates around and one in which it circu- 
lates. The motion starts as libration around with 
mean eccentricity 0.08 and small amplitude (small 
loop on the right) but the system encounters a sad- 
dle and can librate around with large amplitude 
changes of the eccentricity. The full motion is li- 
bration of Am around but it encompasses two 
separate regions. However, every once in a while 
the separatrix is crossed and Am also circulates. 
This suggests the system is on the edge of a reso- 
nance. The motion for planet c has also changed: 
previously it showed a single loop but now it en- 
compasses a moon-shaped regi on that is typical for 
libration in a resonance (e.g. iMurrav &: Dermott 
19991 ) . In the lowest panels Am circulates because 



the planets have crossed the separatrix. The apsi- 
dal no longer prevents them from encounters and 
the system is unstable. 

We searched our solutions for libration of the 
angles a but we found no case where these angles 
librated. We checked this by calculating the reso- 
nant angles for each simulation and tabulating the 
mean, minimum and maximum values. If libra- 
tion occurs then the minima and maxima should 
be different from and 360 and the mean should 
be zero or 180. However, in each case we always 
found the nnnima and maxima to be and 360, 
and the mean to be close to 180, suggesting that 
the angles circulate. We have inspected a large 
sample by eye and conclude that we did not wit- 
ness any libration occurring. While this procedure 
should break down for unstable cases, configura- 
tions which are stable for the entire simulation du- 
ration should have their minima and maxima con- 
fined to a small range. We believe this could be 
caused by the infiuence of GR, which primarily af- 
fects the pericentre precession of close-in planets. 
In turn, this can help break or support a resonant 
configuration and al so affect the stability of multi- 
planet systems (e.g. IVeras fc Fordll2010h . 

The stability dependence of the system on the 
eccentricities of both planets is displayed in Figs.|9] 
and [TOl For these simulations we set the initial 
value of Am = 10° so that we would not witness 
any instability related to this quantity. We var- 
ied ttb and Cb for Fig. [5] and a, and e, for Fig. 1101 
As one may see, the system is mostly stable as a 
function of Pc/Pb and Cb unless Cb ^ 0.2, apart 
from some regions at low period ratio. For such 



14 






pm^ 



A- ■- _. 



C.I ■- .--:■■ -. ■ ■ . ■■■'■W 



■> ^-yi 








y; 



1.3955 1.396 1.3965 1.397 1.3975 1.398 1.3985 1.399 1.3995 



Fig. 9. — Stability map of the HD41248 system as 
a function of the period ratio (horizontal axis) and 
Cb (vertical axis) . The scale on the right indicates 
the value of the stability index D. Blue regions 
are stable, red are unstable. We varied at and e;,. 
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Fig. 10. — Same as Fig. \TU\ but now we varied Qc 
and Be- 



high values of Cf, the apsidal alignment does not 
protect the planets from encounters because the 
orbit of planet b already crosses that of planet c, 
regardless of the eccentricity of planet c. There is 
some fine structure visible in the figure of which 
the line going from (1.3955,0.08) to (1.4,0.03) is 
the most interesting. Below this line there is a 
large region of stable motion and here Aw circu- 
lates. However the eccentricities are small enough 
that the planets do not cross their orbits and they 
are protected by the 7:5 resonance. Above this 
region the stability is guaranteed by the apsidal 
alignment (apart from at high eccentricity). 

In Fig. [10] we see a large unstable region at low 
eccentricity and low period ratio. Here the planets 
are at the edge of the resonance which results in 
chaotic motion and large excursions in eccentric- 
ity resulting in encounters. The stable region on 
the right has the motion being dominated by small 
eccentricity and by the resonance so there are few 
encounters, while at high eccentricity and low pe- 
riod ratio the planets are not in resonance but the 
apsidal alignment prevents close encounters. 

In summary, a system of two super-Earth plan- 
ets is most likely in the HD41248 data when one 
considers the nominal orbits, provided that the 
apses are nearly aligned. Given that Wc appears 
unconstrained by observations, this configuration 
is a viable outcome. The apsidal alignment pro- 
vides protection against close encounters at high 
and low eccentricities, as long as the eccentrici- 
ties are of comparable value. If the eccentricity 
of one planet exceeds 0.2 while the other is nearly 
zero, the system becomes unstable. However, such 
extreme configurations can be ruled out by the ob- 
servational uncertainties. 

7. Discussion 

The population of low-mass exoplanets tends 
to differ significantly from that of more massive 
planets. Observational evidence for a lack of low- 
mass planets orbiti ng the most metal-rich stars 
([Jenkins et al.l 120131 ) points to fundamentally dif- 
ferent evolutionary properties between low and 
high mass planets. As discussed in the introduc- 
tion, there also appears to be a higher fraction of 
dynamically packed low-mass planetary systems 
with planets close to MMR's. 

Here we report the discovery of a possible 7:5 
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MMR planetary system, consisting of two super- 
Earths orbiting a metal-poor star. Alone this dis- 
covery adds to the number of metal-poor stars 
with low-mass planets and lends more w eight to 
the hypothesis presented in I Jenkins et al.l It could 
also be the first confirmed super-Earth system of 
planets in such a second-order resonant configura- 
tion since libration of the longitude of periastrons 
yield a highly stabilized system, and our data set 
and simulations leave the possibility for libration 
of the resonant angles, and hence resonance. Also 
the Doppler data period ratio differs from the true 
7:5 period ratio by less than I part in 300, a strong 
argument in favour of the MMR. 

Around l/3rd of multi-planet systems discov- 
ered by radi al velocities are fou nd to be close 
to a MMR (JLissauer et al.l l201l[ ) but the frac- 
tion actually in the resonance is significantly 
lower. However, iLissauer et al.l show that the 



number of known systems in a MMR is signif- 
icantly more than that expected from a ran- 
domly drawn population. Predictions from con- 
vergent migration show that around 1% of reso- 
nant systems should last for around a disk life- 
time ( Adams et al.l |2008[) whereas planet-planet 
scattering me chanisms yield MMRs around 5-10% 
of the time ( Ravmond et al.l |2008[) . If we take 
the Kepler resonant numbers from iLissauer et al.l 
at face value, since the number of false positives 
from multi-transiting systems is expected to be 
signific antly lower than tha t from single transiting 
events ( Latham et al.ll2011 ). then we might expect 
that the IID41248 planetary system was formed 
through the planet-planet scattering mechanism. 
However, it may not be so straight-forward. 



Terquem fc Papaloizoul (|2007l ) studied the mi 



gration of cores through a proto-planetary disk 
that includes mutual interactions between the 
cores and found that the presence of MMRs is 
common. They found that simulations that be- 
gan with a larger outer planetary distribution ra- 
dius, or that started with significantly lower core 
masses, would tend to produce the 7:5 MMR. 
This indicates that such a second-order MMR 
that forms in this way requires longer timescales, 
since reducing the core masses or increasing the 
planetary core distance from the star gives rise 
to a longer evolutionary time, and thus if the 
HD41248 planets were formed in this way we 
might expect they started with low initial core 



masses, possibly by forming late in the disk evo- 
lution, or they started life much further out in 
the disk than their current semi-major axes, or 
most likely a combination of both low core masses 
and long convergent migration. It is interesting 
to note that even though planet-planet scatter- 
ing si mulations produce mo re MMRs, the results 



from iRavmond et al 



( 20081 ) never resulted in a 
7:5 MMR commensurability. 

Finally, we should note the mass ratio be- 
tween these planets. HD41248 b is found to 
be 12.3 M^and planet c is 8.6 M^, giving rise 
to a mass ratio of 0.7. The simulations of 



Terquem fc Papaloizoul (J2007I ) ended with much 
lower mass ratios for the planets that were found 
in the 7:5 configuration. They found mass ratios of 
0.1 and 0.2, however the more massive planet was 
the one closer to the star, similar to the system we 
present here. We do note that finding MMRs that 
have much lower mass ratios is difficult due to the 
inherent difficulties of first discovering very low- 
amplitude signals in radial velocity datasets, and 
also the difficulty of finding low-amplitude multi- 
planet signals when the super-position of the Kep- 
lerians has a strong dominating signal and a much 
weaker and longer period signal. With this in 
mind, we expect more of these systems to emerge 
from the radial velocity data with the addition of 
more Doppler data and better analysis techniques. 
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Fig. 8. — Evolution of e;, — Atu in the left column and Cc — Atu in the right column. The top panels show 
the evolution before the instability, the middle panels are at the stability edge and the bottom panels show 
the evolution for an unstable system. 
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